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ABSTRACT 

The dynamics of two counter-streaming electron-positron-ion unmagnetized plasma 
shells with zero net charge is analyzed in the context of magnetic field generation in 
GRB internal shocks due to the Weibel instability. The effects of large thermal motion 
of plasma particles, arbitrary mixture of plasma species and space charge effects are 
taken into account. We show that, although thermal effects slow down the instability, 
baryon loading leads to a non-negligible growth rate even for large temperatures and 
different shell velocities, thus guaranteeing the robustness and the occurrence of the 
Weibel instability for a wide range of scenarios. 
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1 INTRODUCTION 

One of the unsolved problems in the astrophysics of extreme 
events is the magnetic field generation in gamma-ray bursts 
(GRBs). It is almost universally accepted the hypothesis 
that the internal and the external shocks in relativistic ejecta 
represent the most likely mechanism for 7-ray production, 
with sync hrotron emission a s the main radiation process 
in GRBs iSari fc Piran|[l99^ . thus indicating the presence 
of strong particle acceleration and near-equipartition mag- 
netic fields. Although the latter has been confirmed by syn- 
chrotron radiation measurements, the origin of the required 
magnetic field is still under strong debate (.Co burn fc Bo ggs, 
l2003l:IWaxmaJ2003D . The Weibel instability (IWeibeMlQSg^ T 
which has received an increasing interest in recent years, 
seems t o be a good candidate to explain magnetic field s in 
GRBs f Medvedev fc LoebllToogl : iMedvedev et al.l 12005*). as 
confirmed by three dimensional particle-in-cell simulations 
JSilva et al J 120031: iFrederiksen et alJl20o4 iNishikawa et alJ 
l2003^ . 

In this paper, we examine the dynamics of two counter- 
streaming electron-positron-ion unmagnetized plasma shells 
with zero net charge, scenario that applies to the internal 
shocks of GRBs. Using relativistic kinetic theory we focus 
our discussion on the growth rate (F) of the electromag- 
netic beam-plasma instability or Weibel instability, since it 
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provides an estimate of the maximum magnetic field gener- 
ated by the instability. Indeed, the saturation level bsat of 
the magnetic field is related to F according to bsat oc 
jMedvedev fc LoeblHooil: iDavidson et~al]ll972l : lYang et alJ 
[^9^) ~A similar procedure was used to examine the mag- 
netic field gener ation due the Weibel instab ility in GRB 
external shocks dWiersma &c Achterbersll2004) . In previous 
analytical works on the Weibel instability in GRBs, the rela- 
tivistic plasma is assumed cold and only purely electromag- 
netic transverse modes are excited; moreover, no arbitrary 
mixture of plasma components has been considered. In this 
work, we employ relativistic kinetic theory to include ther- 
mal efi^ects, with arbitrary temperatures, space charge effects 
and an arbitrary mixture of leptons and baryons. 



It is well known that thermal effects can strongly sup- 
pre ss and even shutd own the current filamentation instabil- 
ity iSi^^^^lJ|200 ^ . Fur thermore, a non relativistic anal- 
ysis (|^ou^^^^L||200^ has shown that the component 
of the electric field E along the wave vector k can also be 
important for two inter-penetrating current beams having 
different temperatures. The different rate at which these 
current beams pinch leads to density filaments, thus cre- 
ating a space charge which reduces the growth rate and 
causes ions to respond. It is important to consider such 
effects because the relativistic shells ejected by the cen- 
tral engine have not only a spread in the bulk velocities 
but should also show the presence of an arbitrary mixture 
of electrons, positrons and protons, with varying tempera- 
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ture from shell to shell |R ees fc Meszaro sll992l : lLazzati el alJ 
ll999l:lPanaitescu. Soada fc Meszarojll999h . 

We have studied several scenarios by changing the prop- 
erties of the colliding shells (e.g. particle number density, 
temperature and bulk velocity of the shell) in order to have 
a more complete analysis of the phenomena that arise in 
these scenarios. We show that, although thermal effects slow 
down the instability, the thermal spread does not shutdown 
the Weibel instability. The presence of ions becomes crucial 
to sustain the instability when shells with different veloc- 
ities are considered. It is shown that baryon loading leads 
to a non-negligible growth rate by lowering the threshold 
predicted for pu rely electro magnetic waves in an ele ctron- 
positron plasma iSilva et "81^2002 : Fonseca el al.l2003l) . thus 
guaranteeing the robustness and the occurrence of the insta- 
bility even for colliding shells with large temperatures and 
different mixtures of plasma components. The threshold for 
the occurrence of the Weibel instability, including relativis- 
tic and space charge effects, is also discussed. 

This paper is organized as follows. The theoretical 
model employed in our analysis is described in Section 2. 
In Section 3 we present a discussion of the numerical results 
regarding the growth rate of the Weibel instability and the 
threshold for its occurrence. Finally, Section 4 summarizes 
our main conclusions. 



2 THEORETICAL MODEL 

To motivate our discussion, let us consider a simple model 
for two counter-streaming plasma shell s eject ed by the cen- 
tral engine of the burst (Lazzati cl al] |l999ll . with respect 
to the c entre of mass ( CM) reference frame of the two-shell 
system llMedvedev fc Loeb..l999.'l . Each shell is characterized 
by an arbitrary mixture of electrons, positrons and protons, 
a bulk Lorentz factor 7 and a bulk thermal Lorentz factor 
Itho = (1 - Ptho)~^'^^ being Ptho the thermal velocity. The 
total charge in each shell and the momentum in the CM 
reference frame are always conserved. 

In order to obtain the dispersion relation for waves 
propagating in the z direction with wave vector k= ke^, 
we use the relativistic Vlasov equation a nd Maxwell's equa - 
tions, following the technique outlined bv lSilva et alJ (j2003), 
yielding: 

(tJ^ — k^C^ + Cxxz) [{'Jj'^ — k^C^ + Cyyz) [i^"^ "f -Dz ) ~ DyCzyZ 

CxyZ ^CyxZ (^ ~f~ DyC zxz\ 
~t~ Dx \CyxzC ZyZ ('^ k C -\- Cyyz^ Czxz\ ~ 



with 



Clmn = ^ UJpOjVnj j dp 



UJVi 



LO — kVz 



Di = Yl Lolojm-i I dp 



LUVl 



u) — kvz 



dp, > Foj, 



(1) 



(2) 



(3) 



where the j component of the plasma is described by 
the corresponding unperturbed normalized distribution 



function Foj, the unperturbed density noj, the mass rrij 
and the plasma frequency ujpoj- For purely electromagnetic 
waves propagating with the same wave vector k= ke^ 
(Silva et al...2002,') . only Cxxz, Cyyz, Cxyz, Cyxz appear in 
the linear dispersion relation; the other terms in Q are 
due to the inclusion of space charge effects, i.e. k-E^ 0. 
Equations Q, Q and ^ provide a starting point for 
the analysis of the Weibel instability in the relativistic 
regime including kinetic effects. Calculation of Q and 
^ for a general distribution function is, in general, quite 
difficult and only possible in closed analytical form for a 
reduced number of distribution functions. We choose a 
waterbag distribution, given by foj = najFoj, with Fqj — 

[1/ {2pzOj)] 5 {Px - PxOj) S [py) {e {pz + PzOj) -0{Pz- PzOj)}, 

where pxOj is the momentum in the x direction and pzOj 
is the momentum thermal spread in the z direction of 
the j plasma component, while Q{x) is the Heaviside step 
function. With this choice of the distribution function, 
the dispersion relation describes the Weibel instability in 
the relativistic war m fluid limit, for con ditions far from 
the kinetic regime jTzoufras et alj l2006ll . Such distribu- 
tion function describes a cold plasma shell in the x and 
y directions, propagating along the x direction, with a 
thermal spread along the z direction. Due to the form of 
the distribution function, the Pm-, Pn- and pz-integrations 
required in Q can be carried out in closed analytical form. 
We use the subscripts "1" ("2") for the shell moving in 
the -t-Cii-direction (-ea;-direction) and we assume, without 
loss of generality, that the shell-1 is faster than the shell-2; 
charge neutrality for each shell is also assumed. 

After normalizing time to l/WpOe-,i (where w^q^- j = 
47rnj.-Q le^/rrig- ) and space to the collisionless skin depth 
c/cjpOe-,ii ^nd performing the integration over the distribu- 
tion function (note that Cxyz = Cyxz ~ Czyz = Dy = and 
Czxz = Dx), the dispersion relation Q can be simplifled to 



{UJ^ - fc" + [Cxxz]) {uj^ + [Dz]) - [Dx]^ = 



(4) 



where the terms in square brackets correspond to the coef- 
ficients from and now given by: 
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I0],s Lo"^ - k'^Pl 



thOj,s 



(5a) 



(5b) 



(5c) 



where the sum in s is over the two shells and the sum in j 
is over the three plasma components, and where /3oj is the 
velocity along the x direction and PthOj is the perpendicular 
thermal velocity of the j plasma component, with the stan- 
dard definitions "foj = (l - Plj - I3th0j)~^^^, ^oj = lOjPo], 
complemented by the definitions [l/'yj^ — J dp(Foj/7), 
and with 



2,3 

E 



70j, 



(6) 
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In order to obtain the growth rate of the instabihty, one has 
to find the roots of the dispersion relation ||1J and consider 
the largest positive imaginary part. Due to its complexity, 
a fully analytical solution leads to cumbersome expressions; 
therefore I0J will be solved numerically. 

However, it is possible to study the threshold condition 
for the Weibel instability, since near the threshold \uP\ <C 1. 
The dispersion relation presents unstable solutions whenever 



< 0, which, from (0J, © and (O yields 



(7) 



where 



s,i=i,i 



PL 



7oj, 



hOj,: 



u 



2,3 

E 

s,j=i,i 



E 



70 j, 



thOj,: 



7oj,s Pthoj,i 



(8b) 



(8c) 



thus Kenera,lizinK the threshold condition derived by 
ISilva et alJ i2002l) . The maximum unstable wavenumber, 
kmax, is then given, from ||7J by 



k^ 



+ 



(j-L-g + T) 



{n - g + j'f - A(T^ - H{g - r)) 



(9) 



which corresponds to a real solution whenever 

e-^>^ (10) 

thus establishing the condition for the existence of unstable 
Weibel modes. It is important to point out that the last 
terms on the right-hand side of Q and HUH represent the 
contribution of the space charge effects. 

Space charge effects arise from the filament pinching due 
to the Weibel instability, which establishes a charge imbal- 
ance in the region surrounding each filament, thus leading to 
an electric field parallel to the wave vector k (i.e. k-E^ 0). 
This (first order) longitudinal electric field component tries 
to prevent pinching, making more difficult the occurrence 
of the Weibel instability. Without the space charge term, 
the threshold is equivalent to that obtai ned by assuming a 
purely transverse electromagnetic mode ijSilva et al.ll2nn3l . 
Since the space charge term is non-negative, it slows down 
the growth of the Weibel instability, and it raises the insta- 
bility threshold. Space charge effects are cancelled out by 
shells having the same density, the same velocity and the 
same temperature. 

To examine the different contribution to the space 
charge effects of the species we are taking into account in 
our model, some straightforward calculations can be done. 
In the limit of an almost cold plasma (uthOj ^ 1), rela- 
tivistic shells {uoj >• 1) and using f3thOe+ ~ VMPthOi, it is 
easy to verify by comparison that, for a scenario with two 
colliding shells having the same density, the space charge 
terms for ions and positrons are of the same order, i.e. 



(if /Hi) / /H^+) ~ 1: baryons and positrons affect 
in the same way the threshold of the Weibel instability 
in what solely concern space charge effects. For scenar- 
ios of our interest, the space charge effect term for a cold 
plasma is cancelled fo r shells moving with the same velocity 
llDavidson el alJl2004l . 

If we neglect space charge effects and keep the previous 
approximations (almost cold and relativistic plasma), from 
nun we recover the condition in :Silva et al.. (,2002), g — T > 



0, which with the definition A, 



written as 



r 1 1 

L "Ij J 703 



3^ — 



2,3 

E 



Aj,, > 0. 



(11) 



The contribution Aj of each species is easily identified. Com- 
parison of the contribution to llll l of baryons and positrons 
at the same temperature yields 



A. 



(12) 



thOe + 



Since (ithoi < I3th0e+ (AhOe+ ~ v^/?thOi) the ratio GUl is 
greater than 1: baryons are more "efficient" than positrons 
to guarantee the occurrence and the robustness of the Weibel 
instability when the temperature increases. What is impor- 
tant for the occurrence of the instability is not the species 
temperature, but the species thermal velocity PthOj instead. 
We will explore this aspect in the next section. 



3 DISCUSSION 

To assure charge neutrality in each shell, we consider an 
electron density noe-,si ^ positron density nQe+ ^ with 
Qs = Koe+.s/'T-Oe- ,sj ^^d the conscqueut ion density noi,s = 
(1 — cts) noe- s', the subscript s denotes the reference shell. 
Our calculations are performed in the CM reference frame, 
which naturally imposes the condition that the total lin- 
ear momentum of the system of two shells in this frame 
is zero. The momentum conservation along the x direc- 
tion yields |?;o,i| (Zi?=i 7ojmoj) ^ = |uo,2| (ELi 7oi?nofc) 
where j,k — 1,2,3 stand for electrons, positrons and ions 
respectively. Since all particle species in each shell have the 
same mean velocity (i.e. the velocity of their own shell), 
the latter is independent of the species considered. The 
momentum conservation establishes limiting values for the 
features of the plasma shells in the CM reference frame, 
and for the parameters describing the shells. We define 
the ratio of the absolute values of the shell velocities as 

(Ef=l 70j"I0j)^ _ I'D, 2 
(Efe = i 70fc»"0fc)2 



V 



Since we assume that the 



shell-1 is faster than the shell-2, the value of V can vary 
within the range ^ ^ 1. Denoting M = rrii/me the ion 
and electron mass ratio and = Joe- s^yoi,sM , it is then 
straightforward to obtain an expression for the ratio of the 
positron and electron densities in the shell-2, which reads 

1 Aj" + niA7 _ 
V 



can only vary in the range: 



Since ^ Q2 ^ 1, the value of oi 



Ar 



Ar' 



(13) 
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Figure 1. Dependence of the maximum growth rate of the Weibel 
instabiUty on the electron thermal Lorentz factor for shells having 
same velocity {p^o = 2), same temperature and same electron 
density, for different mixtures of plasma components: ai 2 = 1 
(dashed line), 01^2 = 0.5 (solid gray line), 01^2 = (solid black 
line). 
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Figure 2. Unstable area (p^^o = 2) for different scenarios. Solid 
lines represent the threshold from IIUI : dashed lines represent 
the limits from momentum conservation [see llijl l . The gray area 
indicates the unstable area for V = 0.75 and rj = 1. 



We observe that, according to limitations in (1131 . ai can 
assume all values between and 1 only when the shells have 
the same velocity (i.e. = 1) and the same temperature 
(i.e. = when y = 1), or equivalently when qi =02. 

As typical para meters for the collid i ng sh ells in GRB 
scenar ios, we follow iMedvedev fc Loebl (Il999l) and iPirarJ 
il999l) . who suggest a bulk Lorentz factors yi„t ~ 4, and 
initial electron thermal Lorentz factors 7i^og ~ 2. Since 

-ythOj = {1 - l3thoj)~^^^ and 7oj = (l - /^oj - 
the perpendicular temperature of the shells is then ob- 
tained from the energy-momentum tensor as T± = 



J dpp 



which for our choice of Foj is 



70 J, s 
2 



1 + 



■In 



-j) (|Silva et all l2002l ). We 



always consider that all particles contained in the same 
shell have the same temperature, whereas different shells can 
have distinct temperatures; the shell electron density ratio 
is indicated by rj = ne-oi/'i-e-02- Three different mixtures 
of plasma components are examined: i) electron-positron 
plasma, which means that ai = 1, ii) electron-positron-ion 
plasma, with ai = 0.5, iii) electron-ion plasma, such that 
Qi = 0. 

When the colliding shells have the same velocity (i.e. 
V = 1, with Pj;o = 2 such that yi„t ~ 4), the same elec- 
tron density (i.e. 77 = 1) and the same temperature, the 
space charge effects are cancelled due to the form of the 
space charge term [see lUnH . The depend ence of the maxi- 
mum growth rate of the Weibel instability on the electron 
thermal Lorentz factor is shown in Fig. Q for this scenario, 
^ ai = Q?2 ^ L As expected, the maximum growth rate 
decreases as temperature increases, until the Weibel insta- 
bility shutdowns for all plasma component mixtures. The 
ion presence is fundamental to maintain the occurrence of 
the instability for larger temperatures than the shutdown 
temperature for electron-positron plasma (corresponding to 
7thOe ~ 1.75). As the baryon component increases in the 
plasma mixture, the shutdown temperature for the Weibel 



instability to occur also increases. This result confirms the 
estimate of the previous section about role of baryons to 
sustain unstable scenarios for higher temperatures [cf. (Il2l l1. 

A mean to analyze a large number of the possible sce- 
narios for the occurrence of the Weibel instability for two 
counter-steaming shells having same temperature, depend- 
ing on the parameters V and rj, is provided by Fig.|5| Here, 
the threshold for the occurrence of the Weibel instability de- 
scribed by HlOjl is represented by solid lines; the points below 
of these lines correspond to unstable scenarios. The dashed 
line represents the lower limit established by 11311 (constraint 
associated with CM reference frame). The negative slope of 
these curves is a clear signature of the role played by baryons 
(ai < 1), which allow the occurrence of the instability for 
larger temperatures than those for mixtures of leptons only 
(ai = 1). 

As for the analysis of shells having same temperature 
but different velocity and electron density, we consider the 
shell-1, moving with p^o = 2 (=> 7i„t ~ 4), to be faster 
than the shell-2 (i.e. < 1). Among the scenarios we can 
examine by varying V and ry, let us study a case with shells 
having the same electron density, such that rj — 1. Since 
V < 1, space charge effects are automatically taken into ac- 
count because of (|Hll. As the previous situation without space 
charge effects, a temperature exists above which the Weibel 
instability in the electron-positron plasma cannot develop 
anymore, whereas when ions are included, they lead to a 
small but non-negligible growth rate. This is shown in Fig. 
13 where the growth rate for all wavenumbers is plotted as a 
function of the electron thermal Lorentz factor for V = 0.75 
and ri — 1 kept fixed. The corresponding area of unstable 
points for this set of parameters is the gray filled area in Fig. 
|5| In order to guarantee momentum conservation, we have 
chosen three values of ai — 1, 0.75, 0.5, considering that for 
a = 0.5 the smallest electron thermal Lorentz factor must 
be 7th0e ~ 1.13. The role played by baryons is illustrated 
more clearly in Fig. 0] where the maximum growth rate of 
the Weibel instability is plotted against the electron thermal 
Lorentz factor (the light gray area is the forbidden area to 




Figure 3. Growth rate of the Weibel instability for all wavenum- 
bers as function of the electron thermal Lorentz factor for shells 
having different velocity {pxO = 2, V = 0.75), same electron den- 
sity {rj = 1), same temperature, for different mixtures of plasma 
components; (a) oi = 1, (b) ai = 0.75, (c) ai = 0.5. 

guarantee momentum conservation for a = 0.5). By com- 
paring Fig. and Fig. |1] for the scenarios with ai — 1 and 
oi = 0.5, even though the inclusion of space charge effects 
leads to a less robust instability (i.e., as the temperature in- 
creases, the instability shutdowns earlier) by increasing the 
threshold for its occurrence, the baryon loading is always 
crucial to the occurrence of the instability for large temper- 
atures. 



4 CONCLUSIONS 

In summary, we have analyzed the possible role of the 
Weibel instability in magnetic field generation in GRB 
internal shocks. Starting from relativistic kinetic theory, 
we have studied the dynamic of two counter-streaming 
electron-positron-ion unmagnetized plasma shells with zero 
net charge. We have examined the Weibel instability in a 
hot plasma with arbitrary mixtures of plasma components 
and we have taken into account the space charge effects, 
so the mode is not purely electromagnetic and transverse 
(k-E^ 0). Our discussion has been focused on the growth 
rate F of the Weibel instability, since the maximum value of 
the magnetic field, which corresponds to the saturation level 
bsat, is related to F according to haat oc F^. Our calculations 
are performed in the CM reference frame, which means that 
the total linear momentum of the two shell system in this 
frame is zero. 

Thanks to an analysis of the threshold for the occur- 
rence of the instability for shells at the same temperature, 
it is possible to locate unstable scenarios by varying shell 
velocities and mixture component densities, taking momen- 
tum conservation into account. For shells with same veloc- 
ity, temperature and electron density (i.e. no space charge 
effects), we have found that, even if the thermal spread slows 
down the Weibel instability and eventually shutdowns the 
instability, the presence of ions guarantees the occurrence 
of the instability for larger temperatures than the shut- 
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Figure 4. Dependence of the maximum growth rate of the Weibel 
instability on the electron thermal Lorentz factor for shells hav- 
ing different velocity {j>xO = 2, V = 0.75), same electron den- 
sity (rj = 1), same temperature, for different mixtures of plasma 
components: ai = 1 (dashed line), ai = 0.75 (solid gray line), 
oi = 0.5 (solid black line). The light gray area indicates the 
forbidden area to guarantee the momentum conservation for 
ai = 0.5. 

down temperature for mixtures of leptons only. In the case 
of shells having different velocity, different electron density 
and same temperature, the occurrence of the Weibel insta- 
bility for large temperatures and consequently the magnetic 
field generation, are still guaranteed by the baryon loading, 
even when the space charge effects are included. Baryons 
lead again to a small but non-negligible growth rate over 
the shutdown temperature for which the electron-positron 
plasma cannot allow anymore the occurrence of the Weibel 
instability. Preliminary PIC simulations qualitatively con- 
firm our theoretical results concerning the growth rate and 
the threshold condition, and will be presented elsewhere. 
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